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Abstract — We  consider  the  problem  of  estimating  the  locations 
of  a  set  of  points  in  a  ^-dimensional  euclidean  space  given  a  subset 
of  the  pairwise  distance  measurements  between  the  points.  We 
focus  on  the  case  when  some  fraction  of  these  measurements  can 
be  arbitrarily  corrupted  by  large  additive  noise.  Given  that  the 
problem  is  highly  non-convex,  we  propose  a  simple  semidefinite 
programming  relaxation  that  can  be  efficiently  solved  using 
standard  algorithms.  We  define  a  notion  of  non-contractibility 
and  show  that  the  relaxation  gives  the  exact  point  locations  when 
the  underlying  graph  is  non-contractible.  The  performance  of  the 
algorithm  is  evaluated  on  an  experimental  data  set  obtained  from 
a  network  of  44  nodes  in  an  indoor  environment  and  is  shown 
to  be  robust  to  non-line-of-sight  errors. 

Keywords :  Non-Line-of-Sight  localization,  semi-definite 
programming,  robust  matrix  decomposition. 

I.  Introduction 

The  problem  of  localization  has  applications  in  many 
interesting  areas  such  as  cyber  physical  systems  [1], 
molecular  biology  [2],  nonlinear  dimensionality  reduction 
[3]  etc.  In  all  these  applications,  pairwise  noisy  distance 
measurements  are  obtained  between  subsets  of  points/nodes 
in  some  Euclidean  space  and  it  is  required  to  estimate  the 
locations  of  these  points  with  high  accuracies.  Given  the 
range  of  applications,  the  problem  has  received  considerable 
attention  in  the  last  decade  and  is  an  active  area  of  research. 
Even  though  the  problem  statement  looks  deceptively  simple, 
the  general  case  is  shown  to  be  NP-complete  [4],  i.e.  finding 
a  valid  configuration  of  points  satisfying  a  subset  of  pairwise 
distance  measurements  is  not  solvable  in  polynomial  time. 
Significant  research  work  is  devoted  to  obtaining  convex 
relaxations  for  the  underlying  optimization  problem  and 
obtaining  conditions  under  which  the  problem  can  be  solved 
in  polynomial  time.  Semi-definite  programming  (SDP) 
relaxations  have  been  proposed  [5],  [6],  [7]  and  the  relaxation 
has  been  shown  to  obtain  the  exact  solution  for  certain  classes 
of  graphs  known  as  uniquely  localizable  graphs,  examples  of 
which  include  d-lateration  and  random  geometric  graphs. 

Most  of  the  existing  literature  on  convex  relaxations  focus 
on  the  case  where  the  pairwise  proximity  measurements 
are  either  known  noiselessly  or  known  with  some  slight 
perturbations  (i.e.  line-of-sight  (LOS)  localization).  However 
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(a)  Sensor  Localization. 


(b)  Vehicular  networks. 


Figure  1.  (a)  Sensors  placed  in  a  field  for  monitoring  the  region,  (b)  High 

accuracy  localization  for  intelligent  transportation  systems. 

in  practice,  a  significant  fraction  of  the  data  could  have 
outliers.  For  example,  in  the  case  of  sensor  localization, 
non-line-of-sight  (NLOS)  propagation  introduces  large  errors 
in  the  distance  measurements.  Spin-diffusion  phenomenon 
in  NMR  spectroscopy  for  protein  molecular  conformation  is 
known  to  produce  bias  errors  in  the  measurements  [8].  Data 
in  higher  dimensions  do  not  all  lie  on  the  manifold  and  there 
could  be  outliers  that  lie  far  away  from  the  surface  of  the 
manifold. 

Traditional  methods  that  do  not  account  for  the  large 
bias  errors  can  result  in  significant  localization  errors.  The 
hardness  of  the  problem  arises  from  the  fact  that  the  bias 
errors  are  typically  large  and  it  is  not  known  a  priori  which 
measurements  are  biased.  There  has  been  some  work  in  the 
literature  that  addresses  this  case  using  convex  relaxations 
[9].  However  these  relaxations  assume  that  the  locations  of 
the  errors  are  known  which  makes  the  problem  considerably 
simpler  to  solve.  Our  focus  here  is  to  address  the  problem 
of  localization  in  the  presence  of  large  unknown  bias  errors. 
For  simplicity,  we  will  keep  the  sensor  localization  problem  in 
mind  as  we  progress  through  the  rest  of  the  paper,  though  all 
the  techniques  can  be  easily  adapted  to  other  applications.  We 
propose  a  semi-definite  programming  (SDP)  relaxation  for  the 


case  of  NLOS  localization  and  provide  conditions  under  which 
the  algorithm  retrieves  the  node  locations.  The  performance 
of  the  algorithm  is  evaluated  through  simulations  and  also 
validated  on  a  data  set  collected  through  a  real  world  indoor 
experimental  setup  of  44  nodes  that  is  publicly  available  [10]. 
The  proposed  relaxation  is  shown  to  mitigate  the  effect  of  the 
NLOS  errors  to  a  large  extent  by  obtaining  accurate  estimates 
of  the  node  locations. 

II.  Problem  Setup 

Consider  a  static  random  placement  of  N  nodes  in  a 
^-dimensional  euclidean  space.  The  locations  of  these  nodes 
are  unknown  and  these  nodes  would  hence  be  referred  to  as 
agents.  We  also  have  M  special  nodes  whose  locations  are 
exactly  known  called  as  anchors.  Nodes  within  a  radius  r  of 
each  other  obtain  relative  distance  measurements  that  could  be 
corrupted  by  noise.  Based  on  these  distance  measurements, 
the  problem  is  that  of  determining  the  locations  of  the 
agents.  Each  of  these  measurements  is  modeled  as  either 
a  LOS-dominated  signal  or  a  NLOS-dominated  signal  by 
choosing  the  observation  noise  to  be  drawn  from  a  mixture 
of  two  distributions. 

Let  Xi  G  Rk  denote  the  location  of  the  ith  agent.  Let  dij  = 
|| Xi  —  Xj\\  denote  the  actual  distance  between  the  ith  and  jth 
node,  denotes  the  corresponding  distance  measurement.  We 
will  model  the  distance  measurement  as  d2j  =  dt2-  +  bij  +  nij , 
where  is  a  zero  mean  additive  noise  (e.g.  gaussian)  and 
b^  is  given  by, 

b^  =  0  if  LOS 
>  0  otherwise. 

Under  NLOS,  b^  is  usually  taken  to  be  a  positive  random  vari¬ 
able  to  model  the  statistics  of  multipath  noise.  The  noise  model 
is  motivated  by  the  time- of- arrival  measurement  modality  for 
distance  measurements.  Since  the  primary  multipath  signal 
travels  a  longer  distance  than  the  true  distance,  the  time  delay 
of  arrival  for  the  NLOS  path  is  larger  than  the  LOS  path  and 
hence  the  NLOS  bias  can  be  modeled  as  additive.  Similarly, 
the  gaussian  noise  component  is  principally  due  to  thermal 
noise  in  the  receiver  and  so  we  take  that  term  to  be  additive 
as  well.  We  will  assume  that  a  fraction  of  the  measurements 
are  NLOS.  We  discuss  a  semi-definite  programming  approach 
to  tackle  this  problem  in  the  sections  to  follow. 

III.  SDP  Relaxation 

Given  the  problem  statement,  we  can  define  a  network 
(A,  A,  £xi  &A-)  D),  of  N  agents  and  M  anchors,  where  X  = 
[xj ;  X2]  •••;  xjf\  G  RNxk  is  the  matrix  of  agent  locations, 
A  =  [aj ;  G  RMxk  is  the  matrix  of  anchor  locations, 

Ex  is  the  set  of  edges  between  agents,  Ea  is  the  set  of  edges 
between  agents  and  anchors  and  D  is  the  matrix  of  pairwise 
distances  between  nodes  that  belong  to  the  edge  set.  An  edge 
is  present  between  two  nodes  in  the  network  iff  we  have  a 
distance  measurement  between  them  (e.g.  if  they  are  within  a 


Algorithm  1  SDP  relaxation  for  NLOS  localization 

l:  Input:  distance  measurements  d^,  (■ i,j )  G  Ex  U  Ea, 
anchor  locations  {a^}^l1  dimension  k. 

2: 


minimize 

x.y.b 


subject  to 


£  ( 


e  iYea  -  d2A  +  b^ 


i,j££x 


£  ([• 


,  T  T I 
e*  aj\ 


Y  X 
XT  Ik 


-  dfj  +  bij 


bij  >  0,  i,  j  £  Ex  U  £a 


Y  X 
XT  Ik 


y  0. 


3:  Output:  minimizer  X. 


distance  r  of  each  other).  We  formulate  the  problem  of  esti¬ 
mating  the  node  locations  as  an  optimization  problem.  For  this 
purpose,  we  will  treat  the  biases  6^  ’s  in  the  measurements,  as 
parameters  that  need  to  be  estimated.  The  maximum  likelihood 
formulation  for  estimating  the  node  locations  as  well  as  the 
biases  gives  us  the  following  optimization  problem, 

maximize  p  ({^}|{xi}£i,  {M) 

subject  to  b^  >0,  i,  j  G  Ex  U  £ a- 


The  above  optimization  reduces  to  the  following  least 
squares  minimization  if  we  take  n^  ’s  to  be  i.i.d  gaussian. 


minimize 

X,B 


E 

i,j<E£x 


\xi  -  xj\\l -dij  +  bij 


V  4 
/ 


E 


\xi  -  aj\\l -d2ij  +  bij 


(2) 


subject  to 


v  ^  J 

bij  >  0,  i,  j  ££xUSa 


where  B  =  {bij}  is  the  matrix  of  bias  variables.  The  above 
minimization  problem  is  highly  non-convex  and  has  multiple 
minima.  We  aim  at  obtaining  a  convex  approximation  to  the 
above  problem. 


The  distance  between  the  agents  can  be  expressed  as  \  \xi~ 
xj\\l  =  eJjXXTeij  where  is  a  vector  in  RN  with  a  1 
and  —1  at  the  ith  and  jth  locations  respectively  and  zeros 
elsewhere.  Using  the  substitution  Y  =  XXT,  we  get  \\xi  — 
Xj  1 1 1  =  efjYeij.  The  minimization  problem  (2)  can  thus  be 


written  as, 


Agent  positions  in  3D  with  the 
same  graph  structure  and  inter 
node  distances 


_  ,  ~  v  2 

minimize  ^  (eJjYeij  -  d2^  +  bij)  + 

’  ’  i,jesx 

_  ,  v  2 

T;  (ef  ye*  —  2ef  Xa^  +  aj dj  -  d j-  +  djj 
i,jesA 

subject  to 

bij  >  0,  i,  j 

y  =  uT, 

(3) 

where  is  a  vector  with  a  one  at  the  ith  position  and  zeros 
elsewhere. 


Since  the  constraint  Y  =  XXT  is  not  convex,  we  relax 
it  to  y  >:  XXT .  Using  Schur  complements,  we  obtain 
the  minimization  problem  as  shown  in  Algorithm  1.  The 
above  semi-definite  program  can  be  efficiently  solved  using 
standard  convex  optimization  techniques.  One  of  the  main 
advantages  of  this  algorithm  is  that  we  do  not  need  to  input 
noise  parameters  to  evaluate  the  optimization  problem,  which 
makes  the  algorithm  useful  in  a  practical  setting. 


Let  us  develop  an  intuition  for  the  relaxation.  Without  the 
relaxation,  in  problem  (3),  we  were  searching  for  a  feasible 
set  of  points  in  the  k -dimensional  space  that  minimize  the 
cost  function.  The  constraint  that  the  points  need  to  be  in 
the  /c-dimensional  space  essentially  makes  the  problem  non- 
convex  and  hard  to  solve.  The  condition  Y  >z  XXT  relaxes 
the  search  for  the  points  in  any  A;' -dimensional  space,  such  that 
k'  <  N.  This  can  be  seen  as  follows.  If  k'  is  the  rank  of  Y, 
then  y  >:  XXT  implies  that  there  exists  X'  G  RNx(k'~k) 
'  XT  " 

X'T 


such  that  y  =  [XX'] 


X  =  [XX']  G  RNxk' 


can  be  thought  of  as  the  set  of  points  in  k'  dimensional 
space  that  are  the  minimizers  of  the  optimization  problem  in 
Algorithm  1 .  Figure  2  illustrates  this  relaxation.  The  price  that 
we  pay  for  this  simplification  is  that,  once  we  get  a  solution 
in  some  k'  dimensional  space,  the  projection  onto  the  original 
fc-dimensional  space  could  yield  bad  results. The  rest  of  the 
paper  explores  conditions  under  which  this  relaxation  actually 
gives  a  solution  in  the  /c-dimensional  space  that  would  be  close 
to  the  true  node  locations. 


A.  Non-cooperative  case 

Consider  the  simple  case  of  a  single  agent  and  multiple 
anchors,  i.e.  N  =  1  and  M  anchors.  The  minimization 
problem  can  be  written  as, 


Figure  2.  Illustration  of  the  SDP  relaxation.  We  search  for  points  in  a  higher 
dimensional  space  satisfying  the  same  distance  measurements  and  the  graph 

tWe  estimated. 

The  following  theorem  states  our  result, 

Theorem  1.  For  the  case  when  the  measurement  noise  is  only 
due  to  bias  errors  i.e.,  rii  =  0,  the  minimization  problem  (4) 
gives  an  exact  solution  as  long  as  the  agent  lies  within  the 
convex  hull  of  two  or  more  anchors  with  which  it  has  LOS 
measurements. 

Proof:  Let  x  be  the  actual  node  location,  and  suppose 
it  is  contained  within  the  convex  hull  of  two  or  more  anchor 
nodes  with  which  it  has  LOS  measurements.  By  contradiction, 
let  x  be  a  solution  to  equation  (4)  such  that  x  7^  x.  Let  W 
be  the  set  of  all  i  such  that  dj  is  a  LOS  measurement  and 
x  is  contained  in  the  convex  hull  of  {a^i  G  W}.  Since  the 
actual  location  of  x  is  a  solution,  the  optimal  objective  value 
is  upper  (and  lower)  bounded  by  0.  Thus  it  must  be  true  that 
\\x  —  ai\\2  <  (dj  =  dj),  Mi  G  W.  For  this  condition  to  hold, 
(x  —  x,cii  —  x)  >  0  ,  Vi  G  W  (consider  3  anchors  with  x  in 
the  middle.  To  be  a  distinct  solution,  x  must  simultaneously 
be  closer  to  all  3  anchors  than  xfs  giving  the  inner  product 
condition).  But  this  is  precisely  the  condition  for  the  existence 
of  a  separating  hyperplane  between  two  convex  sets,  where  the 
sets  are  the  point  x  and  the  convex  hull  of  the  LOS  anchors. 
If  there  exists  a  separating  hyperplane  between  these  two  sets, 
then  x  is  not  in  the  convex  hull  of  {a^,i  G  W},  which  is  a 
contradiction.  So  x  must  be  the  unique  solution.  ■ 

B.  Cooperative  case 

In  this  section  we  consider  the  cooperative  case  where 
we  have  internode  distances  between  the  agents.  Recall  the 
optimization  formulation  in  Algorithm  1.  Note  that  for  every 
term  in  the  cost  function  we  have  a  positive  bias  term  b^, 
that  we  wish  to  estimate  and  use  it  to  offset  any  NLOS  bias 
in  the  corresponding  distance  measurement.  The  optimization 
formulation  does  not  assume  any  prior  knowledge  on  the  bias 
magnitude  or  whether  a  particular  link  is  LOS  or  NLOS. 


minimize 

y,x,{bi} 

subject  to 


lvi  2 

Y(y~  2aix  +  ajai-d%  +  bi'j 

i=  1 

bi  >  0,  V  i 
y  >  xTx 


Consider  a  given  network  of  anchors  and  agents 
=  (X,  A,  £x,£a,  D),  with  an  underlying  graph  of 

connections  between  the  nodes,  where  we  have  an  edge 
(4)  between  two  nodes  iff  we  are  given  a  distance  measurement 
between  those  two  nodes.  Suppose  there  exists  another 
set  of  agent  locations  satisfying  the  same  underlying 


Contracted 


(a)  Contractible  network  (b)  Non-contractible  network 

Figure  3.  (a)  Contractible  networks.  The  dotted  lines  and  corresponding 

nodes  represent  a  contracted  network  with  inter  node  distances  smaller  than 
that  of  the  original  network. (b)  Non  contractible  network.  See  that  the  agent 
locations  cannot  be  shifted  without  increasing  at  least  one  of  the  distances. 


graph  topology  as  the  original  network  but  with  inter¬ 
node  distances  restricted  to  the  edge  set,  shorter  than  the 
original  network  of  nodes  i.e.,  A/2  =  (X,A,£x,£a,D), 
where  Dij  <  D^,V(i,  j)  E  fx  U  £a-  Figure  3(a)  is  an 
example  of  such  a  network.  In  such  a  case,  the  optimization 
problem  cannot  distinguish  whether  the  given  set  of  distance 
measurements  belonged  to  the  original  network,  or  whether 
the  distance  measurements  were  obtained  from  the  second 
network  corrupted  by  positively  biased  noise.  Thus  we  can 
only  hope  to  have  a  unique  solution  as  long  as  the  original 
network  has  a  graph  structure  that  rules  out  the  existence 
of  a  different  set  of  node  locations  satisfying  the  same 
graph  structure  and  smaller  inter-node  distances.  A  fully 
connected  graph  is  one  such  example.  Figure  3(b)  shows 
another  example  of  such  a  network.  Note  that  the  node 
locations  of  the  agents  cannot  be  shifted  without  increasing 
at  least  one  of  the  distances.  Both  these  examples  are  in 
two-dimensions  where  we  assumed  that  the  nodes  are  also 
placed  in  two-dimensions  (i.e.  k  =  2).  However  we  saw  that 
the  SDP  relaxation  allows  us  to  search  in  for  node  locations 
in  a  higher  dimensional  space.  Hence  we  would  need  this 
property  to  hold  in  higher  dimensions  too.  This  motivates  the 
definition  of  non-contractible  networks. 

Definition:  A  network  (X,  A ,  £4,  D)  is  said  to  be  non- 

contractible  if  there  exists  a  unique  X  E  RNxk  satisfying 
the  location  constraints  imposed  by  D  and  there  exists  no 

Xj  E  ffiA,  j  =  1,  2, ..,  N,Wh>  k  such  that 

Wiai\ 0)  —  x'j\\2  <  d2ej  V  (£,j)e£A 

\Wi~x'j\\2  <  d2^  V  ( i,j)e£x 

Xj  /  (xj ;  0)  for  some  j  E  {1,  2,  3, ....,  N} 

Figure  3  shows  examples  of  networks  that  are  contractible 
and  non-contractible  in  M2.  Note  that  even  if  one  node 
lies  outside  the  convex  hull  of  the  anchors,  the  network  is 
contractible.  However  the  converse  does  not  hold. 

The  following  theorem  establishes  that  non-contractibility 


Figure  4.  Error  as  a  function  of  the  fraction  of  NLOS  measurements  under 
the  different  formulations.  The  NLOS  bias  was  uniform  on  [0,6],  the  noise 
was  zero-mean  Gaussian  with  cflos  =  0.02,  and  error  values  are  averaged 
over  10  trials  each. 


is  a  necessary  and  sufficient  condition  for  the  SDP  relaxation 
to  provide  an  exact  and  unique  solution. 

Theorem  2.  For  the  case  when  the  measurement  noise  is  only 
due  to  bias  errors  i.e.  riij  =  0, 

•  Algorithm  1  gives  the  exact  solution  if  the  underlying  LOS 
network  of  nodes  is  non-contractible. 

•  If  the  max  rank  solution  of  the  optimization  problem  in 
Algorithm  1  is  k,  then  the  solution  is  exact. 

Proof:  The  intuition  for  the  proof  relies  on  the  arguments 
made  for  the  definition  of  non-contractibility.  The  sufficiency 
of  the  non-contractibility  condition  is  proved  using  arguments 
of  contradiction.  See  the  appendix  for  a  formal  proof. 

■ 

Remark:  The  notion  of  non-contractibility  is  a  general¬ 
ization  of  the  uniquely  localizable  condition  [5]  for  NLOS 
localization.  However,  the  question  of  what  networks  are  non- 
contractible  in  practice  is  still  an  open  question.  There  has 
been  recent  progress  in  the  LOS  case  in  characterizing  the 
classes  of  networks  that  are  uniquely  localizable  using  graph 
rigidity  theory  [6].  Analogously,  one  could  hope  to  utilize 
results  from  tensegrity  theory  [11],  to  characterize  networks 
that  are  globally  stable  and  non-contractible.  This  is  a  future 
direction  of  research. 

IV.  Validation 

A.  Simulation 

Simulations  were  carried  out  for  the  proposed  SDP 
relaxation  and  compared  against  the  scheme  proposed  by 
Biswas  et.  al.  [5].  We  chose  this  baseline  system  to  compare 
as  we  are  not  aware  of  any  literature  on  convex  relaxations 
that  accounts  for  NLOS  without  assuming  prior  knowledge 
on  the  links  that  are  NLOS.  The  SDP  relaxation  in  [5]  is 
robust  to  small  percentages  of  NLOS  noise  but  fails  for  larger 
percentages  given  that  they  do  not  explicitly  account  for 
NLOS.  The  setup  consisted  of  50  nodes  randomly  placed  in  a 
[—1,1]  x  [—1,1]  grid  with  15  anchors.  Anchors  were  placed 


Figure  5.  Per-node  error  as  a  function  of  communication  radius  on  the  data 
set  obtained  by  Patwari  et  al  [12]. 


Figure  6.  Node  location  estimates  for  the  experimental  data  set  using  our 
SDP  relaxation. 


uniformly  along  the  upper  and  lower  boundaries  of  the  grid. 
The  radius  of  connectivity  was  taken  as  r  =  1.5,  and  the 
gaussian  noise  standard  deviation  was  set  to  0.02.  The  NLOS 
noise  was  taken  to  be  uniform  over  the  interval  [0, 6].  Results 
were  obtained  after  averaging  over  10  trials. 

Figure  4  shows  the  variation  of  the  error  in  the  node  location 
as  a  function  of  the  fraction  of  NLOS  measurements  for  the 
relaxation.  The  error  metric  is  the  per  node  average  mean 
squared  error  in  the  node  location  estimate  i.e.  \\x-x\\/VN. 
It  can  be  seen  that  the  proposed  SDP  relaxation  is  robust  to 
a  large  fraction  of  NLOS  errors. 

B.  Experimental  Results 

The  SDP  relaxation  was  tested  against  data  obtained  from 
real  world  experiments.  We  use  the  data  set  published  by 
Patwari  et.  al.  [12].  The  experiments  were  conducted  in  an 
indoor  office  environment  by  simulating  44  node  locations 
using  a  transmitter  and  a  receiver  and  obtaining  pairwise 
time-of-arrival  measurements.  The  environment  had  a  lot 
of  scatterers  and  nearly  all  the  distance  estimates  have 
strong  NLOS  biases  in  them.  The  authors  in  [12],  artificially 
subtracted  out  the  NLOS  biases  and  obtained  the  node 


Figure  7.  Node  location  estimates  for  the  experimental  data  set  using  the 
scheme  by  Biswas  et  al.  [5]. 

location  estimates  using  LOS  algorithms.  They  report  a 
per  node  localization  error  of  1.26m  after  utilizing  all  the 
measurements  i.e.  a  fully  connected  network.  We  use  the 
same  data  set  without  subtracting  out  the  NLOS  biases  and 
show  that  the  SDP  relaxation  performs  very  well  even  in  this 
setup.  The  node  location  estimates  are  obtained  using  only 
a  subset  of  measurements,  based  on  an  arbitrarily  defined 
radius  of  communication.  The  SDP  relaxation  does  not  need 
any  inputs  regarding  the  noise  parameters  which  is  required 
by  most  of  the  existing  algorithms  in  the  literature. 

Figure  5  shows  the  per  node  error  as  a  function  of  the 
communication  radius.  The  per  node  error  is  seen  to  be  less 
than  a  meter  even  when  the  NLOS  biases  are  included  and 
only  a  subset  of  the  measurements  are  used.  However  the 
performance  declines  as  the  radius  becomes  too  large  or  too 
small.  The  reasoning  is  that  for  small  radii,  each  node  may 
not  have  enough  information  to  fully  resolve  its  position  and 
the  network  could  become  non-contractible.  For  large  radii, 
the  multipath  bias  is  more  likely  to  be  large,  causing  errors 
which  can  be  reasoned  as  follows.  The  measurement  model 
that  we  considered  has  the  bias  added  to  the  square  of  the 
distances.  However  in  practice  the  bias  would  be  added  to  the 
actual  distance  value.  Hence  when  we  square  the  distances 
and  apply  our  measurement  model,  the  bias  values  would  be 
a  function  of  the  distance  and  hence  would  increase  as  the 
distances  increase.  This  partially  explains  the  reason  why  the 
performance  degrades  slightly  as  we  increase  the  radius  of 
connectivity.  Figure  6  and  Figure  7  show  the  estimated  and 
true  node  locations  for  the  experimental  data  set  obtained  from 
the  SDP  relaxations. 

V.  Conclusion 

In  this  paper,  we  considered  the  problem  of  estimating  the 
locations  of  a  set  of  points  in  a  fc-dimensional  euclidean  space 
given  pairwise  distance  measurements  amongst  the  points.  We 
considered  the  case  when  some  fraction  of  the  measurements 
could  be  arbitrarily  corrupted,  proposing  a  SDP  relaxation 
formulation  for  the  problem.  The  relaxation  was  shown  to 
give  the  exact  solution  when  the  underlying  graph  is  non- 
contractible.  Simulations  and  testing  on  real-world  data  shows 


that  the  algorithm  is  quite  robust  to  NLOS  noise.  Future 
work  involves  strengthening  the  theoretical  results  in  order 
to  provide  meaningful  guarantees  in  practice.  It  would  be 
interesting  to  characterize  the  classes  of  graphs  that  are  non- 
contractible.  We  also  need  to  theoretically  prove  that  the  SDP 
relaxation  is  robust  to  thermal  noise. 
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solution  X ,  there  exists  a  solution  X'  of  rank  higher  than  k 
that  satisfies  Y  >z  X,T X' .  X'  satisfies 

Y-XTX  =  XfTXf. 


Define 


X 


We  then  have 

XTX 


X 

X' 


XTX  +  X'TX’ 

Y. 


Let  Elos  be  the  edge  set  of  the  LOS  links. 

eJjYeij  -  +  bij  =  0, 


eTYe ••  <  d2 

C IJ 1  3  — 


Xi  -  Xi 


<  4 


,  Vi,  j  G  Elos  Li  Ex ? 
2  Vi,  j  G  Elos  n  Ex- 


We  also  have, 

[ef-aj] 


r  y  x] 

1 

* 

~at 

ejYei  -  2aJ Xei  +  ajat  <  d2u, 


eJXTXei  -2[aj0] 


X 

X' 


2 

ll(®*;0)  -  LW2 

for  Vi,  I 


< 

4 

< 

d2 

aio 

< 

d2 

aiO 

< 

4 

Elos  Li  £a 


Thus,  we  just  showed  that  the  underlying  network  is  con¬ 
tractible  which  is  a  contradiction  to  our  initial  assumption  that 
the  network  is  non-contractible. 


Appendix  A 
Proof  of  Theorem  2 

Lemma  1.  If  the  max  rank  solution  of  the  optimization 
problem  in  Algorithm  1  is  k,  then  it  is  unique. 

Proof:  We  will  prove  this  by  contradiction.  Suppose 
(Xi,Yi,5i)  and  (X2,  F2,  B2)  are  two  solutions,  then  for 
any  /?  <  1  one  can  easily  check  that  /?(Xi,  Yi,  Bf)  +  (1  — 
/3)(X2,  Y2,B2)  is  also  a  solution.  Since  the  max  rank  is  k,  we 
should  have  that, 

m  +  (1  -  f3)Y2  1 3X T  +  (1  -  p)Xl  ‘ 

_  I3X1  +  (1  -  P)X2  4 

be  rank  k.  This  gives  us 

pY1  +  (l-p)Y2  =  (pXi  +  (1  —  (3)X2)t(/3Xi  +  (1  —  fi)X2). 

where  Y\  =  X^Xi  and  Y2  —  X^X2.  Rearranging  terms,  we 
get  \\Xi  —  X2||2  =  0  which  implies  X\  =  X2.  ■ 

We  will  now  show  that  if  the  underlying  LOS  network 
is  non-contractible  then  the  solution  is  unique.  Assume  that 
each  node  has  at  least  one  LOS  measurement.  We  will  prove 
uniqueness  by  contradiction.  Suppose  that  the  network  is 
non-contractible  and  the  solution  is  not  exact.  Then  by  the 
definition  of  non-contractibility  and  by  the  non-uniqueness  of 


